SU(4) lattice gauge theory with decuplet fermions: Schrodinger functional analysis 
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We complete a program of study of SU(A'^) gauge theories coupled to two flavors of fermions in 
the two-index symmetric representation by performing numerical simulations in SU(4). The beta 
function, defined and calculated via the Schrodinger functional, runs more slowly than the two-loop 
perturbative result. The mass anomalous dimension levels off in strong coupling at a value of about 
0.45, rendering this theory unsuitable for walking technicolor. A large-A'^ comparison of this data 
with results from SU(2) and SU(3) reveals striking regularities. 
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I. INTRODUCTION 

In the last few years there has been an explosion of in- 
terest in lattice simulations of theories with gauge fields 
coupled to a large number of fermionic degrees of free- 
dom, either many flavors of fermions in the fundamen- 
tal representation or a few flavors of fermions in higher- 
dimensional representations of the gauge group [H- The 
primary motivation to study these models is their po- 
tential use as technicolor theories — extensions to elec- 
troweak theory with no fundamental Higgs field i2] ■ To 
be a viable candidate, a theory must exhibit spontaneous 
symmetry breaking, providing Goldstone bosons to be 
eaten by the electroweak gauge bosons. It must also have 
a slowly running gauge coupling and a mass anomalous 
dimension 7™ near unity, in order to give realistic masses 
to Standard Model fermions while correctly suppressing 
flavor-changing neutral currents. 

Expectations for candidate theories, based on the two- 
loop beta function, are easily stated [1, 0]. Too many 
fermionic degrees of freedom render the Gaussian fixed 
point infrared stable. With a smaller number of fermionic 
fields, asymptotic freedom returns but the beta function 
possesses a zero, signaling the presence of an infrared- 
attractive fixed point (IRFP). Theories with IRFP's 
are said to reside within a "conformal window" in the 
{Nc,Nf) plane, wherein the IR physics displays confor- 
mal invariance and no particle spectrum. The fermion 
mass is then a relevant parameter whose presence in- 
duces a mass gap. Further decreasing the number of 
fermion fields takes us out of the conformal window to 
the QCD-like domain of confinement and chiral symme- 
try breaking. Just outside the conformal window, there 
may be a borderland where the beta function approaches 
zero without actually crossing it; this gives candidates for 
"walking technicolor," where the running coupling comes 
to a near-standstill for many decades in the energy scale, 
until chiral symmetry breaking eventually sets in. 

This paper presents a study of the running gauge cou- 



pling and mass anomalous dimension for SU(4) gauge 
fields coupled to two flavors of decuplet fermions, that 
is, fermions in the two-index symmetric representation. 
It is a continuation of our previous work on SU(2) and 
SU(3) gauge theories with fermions in the corresponding 
representations [Bl-fTlj. Recent interest in these models 
dates from the proposal by the authors of Refs. (l^ . [Tsj 
that they might make good technicolor candidates. All 
three theories possess two-loop IRFP's, raising the pos- 
sibility that they might become walking theories when 
studied nonperturbatively. In addition, the small num- 
ber of flavors was seen as favorable from the point of view 
of evading precision electroweak constraints. 

Using the Schrodinger functional (SF) method [13413, 
we were able to confirm an IRFP in the SU(2) theory, 
placing it within the conformal window @; our most re- 
cent result for the SU(3) theory was the same, but at a 
lower level of statistical confidence [ll| . Our calculations 
of 7m allowed a strong claim, in both cases, that jm lev- 
els off at strong coupling so that it never exceeds 0.45; 
thus neither theory can be used for walking technicolor, 
regardless of the existence of the IRFP. In an effort to 
find a phenomenologically viable theory with 7™ ~ 1, we 
turn to the SU(4) theory. Our motivation for this lies in 
the one-loop expression for 7^, 
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which is proportional to the quadratic Casimir operator 
C2{R) of the fermion representation; in going from the 
sextet SU(3) theory to decuplet SU(4), C2[R) increases 
from 10/3 to 9/2, or 35%. If the non-perturbative result 
were to follow this pattern then a value near 1 would be 
within reach. 

As in our earlier work, we apply the SF method to cal- 
culate the nonperturbative beta function. The original 
method yields a discrete analogue of the beta function 
that gives the change in the running coupling when the 
length scale changes by a fixed ratio. We have noted @ 
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that when the coupHng runs slowly, as it does in the the- 
ories at hand, an approximate result for the usual beta 
function can be obtained directly. SF techniques also al- 
low one to measure 7„i from the volume dependence of 
Zp^ the renormalization factor of the pseudoscalar den- 
sity 'ipT'^lS'lp-^ 

Analyzing these theories via numerical simulations 
presents a different set of problems than is seen in 
QCD [31]. The major new feature is the slow running 
of the coupling. As we have noted, the slow running 
simplifies much of the analysis, both of the beta func- 
tion and of 7m . On the other hand, it is difficult to tell 
a slowly running coupling from one that does not run 
at all, so that the location of any IRFP may be poorly 
determined. The main problem with slow running, how- 
ever, is that if the bare coupling is tuned to be strongly 
interacting at long distance, then it will be strongly in- 
teracting at short distance. This raises the possibility of 
strong lattice artifacts, in particular the appearance of 
unphysical phase transitions. 

Such a transition is present in all our candidate theo- 
ries. Our lattice simulations use fermions with a Wilson- 
type discretization, which breaks chiral symmetry. Simu- 
lations are done at zero quark mass, but the quark mass 
is a derived quantity, determined from the axial Ward 
identity. Reaching zero fermion mass involves tuning the 
hopping parameter at fixed bare coupling. At strong 
coupling a lattice transition occurs, and when it does, 
nowhere does the fermion mass vanish. Instead, it jumps 
discontinuously from positive to negative value. The ab- 
sence of a massless theory in strong coupling makes it 
impossible to apply the SF method there. To evade this 
problem, we change the lattice discretization to push the 
transition away from the region of bare parameter space 
where we wish to run. 

Our result for the beta function is similar to what we 
found in the triplet SU(2) and sextet SU(3) theories. The 
nonperturbative result is consistently smaller in magni- 
tude than the two-loop estimate, but in this case, as in 
SU(3), we cannot state definitely that it crosses zero. 
The result for 7^ follows closely the pattern of the other 
two theories: It follows the perturbative line, Eq. ([TJ, 
for weak coupling but departs from it and saturates at 
7m ~ 0.45 (see below for error estimates). 

The outline of the paper is as follows. In the next 
section we describe our generalized lattice action, formu- 
lated to suppress lattice artifacts in the strong-coupling 
region. In Sec. IIIII we present our results for the beta 
function and in Sec. IIVI the mass anomalous dimension. 
Since we now hold results for the SU(iV) theories with 
iV = 2, 3, and 4, we can discuss them in the language of 
large- gauge theories; we show in Sec.|V]that the consis- 
tency of the three theories is remarkable. The appendices 



present our method of smearing gauge links in SU(4), 
perturbative and nonperturbative tests of our gauge ac- 
tion, and a tabulation of our simulation ensembles. 



II. GENERALIZING THE LATTICE ACTION 

We study the SU(4) gauge theory coupled to two fla- 
vors of dynamical fermions in the symmetric representa- 
tion (the decuplet) of the color gauge group. Our tech- 
niques are mostly identical to our previous work with the 
SU(2) and SU(3) theories We have aheady used 

a generalized gauge action for the SU(3) theory; here we 
give a more thorough discussion. 

We use the Wilson fermion action with added clover 
term [s^- The gauge connections in the fermion ac- 
tion are defined with a differentiable hypercubic (nHYP) 
smearing [s^ [s^l of the fundamental links, from which 
the decuplet gauge connection for the fermion operator 
is constructed.^ We began this project with the usual 
single-plaquette Wilson gauge action. The parameters 
that are inputs to the simulations are the gauge coupling 
l3 = S/f^g and the fermion hopping parameter k, related 
to the bare mass uiq by k = (8-1- 2mQ)~^. The clover 
coefficient is set to its tree-level value of unity. 

The Schrodinger-functional study of the running cou- 
pling is carried out at zero fermion mass, which defines 
the critical hopping parameter Kc(/3)- We define the 
fermion mass rriq and its related critical hopping param- 
eter Kc through the axial Ward identity (AWI), 
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where the axial current Af, 
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"^It^Toir"- /2)ip, the pseu- 
doscalar density — '4'^^{t°' /2)\p, and is a gauge- 
invariant wall source &i t — a. (See [1] for more details.) 
With the simplest lattice action — the single-plaquette 
gauge action and thin-link Wilson fermions — it turns out 
to be impossible to set the fermion mass to zero when 
the gauge coupling is strong. In that regime there is a 
first-order phase boundary at which the AWI mass jumps 
across zero as k is adjusted, never taking the value zero. 
Thus there is no place where the fermions are massless 
in strong coupling: The Kc line simply terminates. This 
has been observed in theories with many Wilson-type 



fermions — many fundamental flavors [35 - 37 , or a few fla- 
vors of higher-representation fermions [8 , 38 , [39| . For our 
purposes, this phase transition prevents the extension of 
the SF calculation to strong coupling. 

The transition is a lattice artifact. Its location can be 
shifted by changing the lattice action. For triplet SU(2) 
with Nf = 2, SF calculations with thin-link fermions ^25l . 



For other ap 
dates, see 123 



)lications of the SF method to technicolor candi- 



^ The extension of this smearing to the gauge group SU(4) is de- 
scribed in Appendix 1X1 
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FIG. 1: Inverse SF coupling versus gauge coupling /3 for sev- 
eral choices of /3io, measured with short runs on a 6* lattice. 
The connected data sets are for (right to left) /3io = —0.5, 
0, 0.5, 1.0, 1.5. For the first three values of /3io, the vertical 
lines mark the appearance of the first-order transition that 
makes Kc disappear for smaller /3. The horizontal dashed line 
near the bottom of the graph marks the location of the Banks- 
Zaks (two-loop) fixed point. The horizontal dotted line marks 
where the one- loop 7m (ff^) is equal to unity. 

[26j were hampered by this transition; with nHYP clover 
fermions we pushed the transition back and exposed the 
IRFP 0. For SU(3) and SU(4) this change of action 
turns out to be insufficient and the first order transition 
remains at fairly weak coupling. 

We find that a modification of the gauge action can 
move the transition away and allow the study of stronger 
couplings. We supplement the original plaquette term 
with an additional plaquette term, constructed with the 
same link as is used in the fermion action — a fat link in 
the higher representation. The action is 

Sg = ^^ReTr{/^(a:)t/,(a; + A);7^(a; + i>){/t(x) 

+ ^ E ReTrT/^(a;)K(a; + fi)V^ix + i^)Vjix), 

(3) 

where (x) is the thin link and (x) the fat link in the 
decuplet representation. We will refer to this action as 
a "soft gauge action."^ We note that the fat plaquette 



^ We offer a perturbative analysis of soft gauge actions in Ap- 
pendix [B] as well as evidence that they do not negate the advan- 
tages of fattening the links in the fermion action. 



FIG. 2: Inverse SF coupling versus inverse bare coupling l/(/o 
as defined in Eq. Q. Plotting symbols are as in Fig. [T] plus 
additional single points for /3io = 2.0 (star), /3io = 3.0 (x), 
and /3io = 3.5 (triangle at the top). The diagonal line has 
unit slope. 



term would appear in a hopping-parameter expansion in 
O(k^); such an induced term has been blamed for the 
first-order phase boundary (38| . so it makes sense to try 
to adjust its strength. 

The advantages of the soft gauge action can be seen in 
Fig. [1] a plot of the inverse SF coupling 1 /g^ versus (3 
for various values of /3io, measured on a fixed lattice size. 
Our goal is to investigate small values of in partic- 
ular the neighborhood of the Banks-Zaks fixed point and 
of the coupling where the one-loop anomalous dimension 
7m (ff^) is equal to unity. These are marked by the hor- 
izontal lines in the figure. The circles denote results for 
Pio — 0; meaning the original thin- link plaquette action. 
The vertical line just to the left of the leftmost circle 
marks the bare coupling where the first-order boundary 
appears and there is no Kc- Thus one cannot investigate 
the region 1/g^ ^ 0.25 with the plaquette action. In- 
creasing /3io pushes the phase transition to smaller /?; as 
it turns out, this allows us to reach smaller values of 1/g^ 
before encountering the transition. For /3io ^ 1 we can 
no longer find the transition in the coupling region we 
studied. 

Having data from many lattice actions allows us to test 
universality, with the results shown in Fig. [2l We plot 
all the lattice data against the perturbative bare lattice 
coupling. 
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TABLE I: Running coupling, Eq. ((6|, evaluated at the bare 
couplings (/3,/3io,Kc) on lattices of size L. The omission of 
the result for L = 16 at /3 = 5.0 is explained in Appendix O 



p 


/3io 




1/5^ 




L — Qa 


L = 8a 


L = 12a 


L = 16a 


10.0 





0.5778(19) 


0.5440(40) 


0.5241(64) 




9.0 





0.4095(19) 


0.3915(37) 


0.3610(42) 


0.3334(72) 


8.5 





0.3199(19) 


0.2979(32) 


0.2780(69) 




7.0 


1 


0.4377(11) 


0.4217(26) 


0.4024(52) 




6.0 


1 


0.3101(27) 


0.3019(32) 


0.2844(64) 




5.0 


1 


0.2157(27) 


0.2103(40) 


0.2002(53) 


* 


4.5 


1 


0.1714(27) 


0.1722(32) 


0.1721(47) 


0.1631(67) 


4.0 


1 


0.1443(24) 


0.1389(26) 


0.1332(35) 


0.1413(46) 


3.5 


1 


0.1070(18) 


0.1106(29) 


0.1087(42) 





[cf. Eq. (jB3|) ]. The points indeed collapse to a com- 
mon curve, especially in weak coupling. In perturba- 
tion theory, the bare and SF couplings are related by 
5^ = 5§ + C'5o + • • • , or 

l/g^^llgl-C+--- . (5) 

The solid diagonal line is plotted to show that the slope 
of the data indeed approaches 1 in weak coupling. 

III. BETA FUNCTION 
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FIG. 3: Running coupling calculated on lattices of 

size L (see Table H]). The crosses are from simulations with 
/3io = 0: top to bottom, /3 = 10.0, 9.0, 8.5. The circles are 
from simulations with /3io = 1: top to bottom, /3 = 7.0, 6.0, 
5.0, 4.5, 4.0, 3.5. The straight lines are linear fits [Eq. (O] 
to each set of points at given (/3,/3io); the slope gives the 
beta function. The dotted line shows the expected slope from 
one-loop running. 



The computation of the running coupling proceeds as 
described in Ref. [1], with adaptations to the SU(4) case. 
We set boundary conditions on the gauge fields as de- 
scribed in ^25], while the fermions obey the usual homo- 
geneous boundary conditions at t = 0,L. The coupling 
emerges from a measurement of the derivative of the ac- 
tion with respect to a parameter rj in the boundary gauge 
field. 



K 



dSg 

dr] 



tr 



dri Dp 



J7=0 



(6) 

The constant K can be calculated from the derivative of 
the classical continuum action with respect to ry, giving 
K = — 127r. For details of the ensembles generated see 
Appendix O 

Our results for the running coupling are listed in Ta- 
bleland plotted in Fig. [21 It is convenient to define the 
beta function f3{u) for u = as 



dil/g^) 2/3(g2) 



= 2u^l3{l/u) (7) 



dlogL g-i 

in terms of the conventional beta function (3{g'^)- As 
discussed in Ref. [9^1, the slow running of the coupling 
justifies extracting the beta function at each Hc) from 
a linear fit of the inverse coupling 

u{L/a) = Co + cilog-^. (8) 
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FIG. 4: Beta function /3(«) as extracted from the linear 
fits ((HJ, plotted as a function of u{L = 8a). The squares 
are from the /3io = data while the circles are from /3io = 1. 
Plotted curves are the one-loop (dotted line) and two-loop 
(dashed line) beta functions. 
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FIG. 5: Comparison of different fit types. Empty symbols as in Fig.|4l Left panel: full symbols derive from linear fits in which 
the L = 6a points are omitted. Right panel: full symbols derive from fits to Eq. pop . in which a log^ term has been added. 
Results are plotted against u{L = 8a). The filled symbols have been slightly displaced to the right. 



ci is an estimate for the beta function /3 at this coupling. 
These fits are shown in Fig. [3l Each fit was done using 
all the available volumes at the given bare parameters. 
Values of the beta function (3{u) obtained from these fits 
are plotted as a function of u{L = 8a) in Fig. S) There 
is some discrepancy between the results for jSio ~ and 
those for /3io = 1; as we shall see, this discrepancy is 
not robust under changes of the fitting procedure. Also 
shown are the one- and two-loop approximations from 
the expansion 



/3(u) = - 



2bi 
T67r2 



262 



(167r2)2 u 



(9) 



where bi = 20/3 and 62 = -260/3. 

The assumption behind the linear fits is that (3 is small 
so that u{L/a) changes very slowly with the volume; this 
behavior is apparent in Fig. [3] Indeed the fits have good 
X^, which justifies our hypothesis. Corrections to the 
simple model (|8]) come from discretization errors, as well 
as from the slight deviation from constancy of the con- 
tinuum beta function over the range of volumes. Dis- 
cretization errors have the form of powers of a/L. We 
have found that such corrections are only loosely con- 
strained in a generalized fit; thus we prefer to estimate 
the uncertainty due to these corrections by redoing the 
linear fits while omitting the smallest lattice, L = 6a. 
The results are shown in the left-hand panel of Fig. [S] 
While the error bars have increased, on the whole the 



results are stable. 

Deviations from constancy of the (continuum) beta 
function give rise to higher powers of logL/a. Adding 
the next-to-leading term, at each bare coupling we fit 



i{L/a) = Co + ci logL/8a + C2(logL/8a)^ 



(10) 



From the definition of the beta function it follows that ci 
continues to provide an estimate for the beta function at 
u — l/g^{L = 8a). The results of these fits are shown in 
the right-hand panel of Fig. [51 This time there is hardly 
any change compared to the linear fits of Eq. ([5]).^ 

While the data trend towards a zero crossing, meaning 
an IRFP, we cannot confirm the existence of this crossing. 
It is possible that the beta function turns away from zero, 
resulting in a walking scenario. If the function does cross 
zero, we can offer an estimate of the crossing point. For 
each fit type, we determine the zero of the beta function 
from a linear fit of /? vs. u, using the /3io = 1 points. 
In fitting to any of Figs. HHSl we find little difference 
whether we fit to 4, 5, or 6 points starting at the lowest 
u. Moreover, the three figures give nearly equal central 
values for the crossing, and comparable error bars. The 
largest la interval covers all the others. Taking it as our 



* Note that dropping L = 6a for the cases where there are only 
three volumes leaves no degrees of freedom for the linear fit. 
Here, too, fitting the cases with only three volumes leaves no 
degrees of freedom for the fit. 
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TABLE II: Pseudoscalar renormalization factor Zp evaluated 
at the couplings (/?, /3io, Kc) for lattice sizes L. 



/3 






Zf 






L = 6a 


L = 8a 


L = 12a 


L — 16a 


10.0 





0.2210(3) 


0.2011(6) 


0.1754(7) 




9.0 





0.1852(4) 


0.1618(4) 


0.1324(4) 


0.1168(6) 


8.5 





0.1565(4) 


0.1331(5) 


0.1066(7) 




7.0 


1 


0.1869(2) 


0.1682(3) 


0.1473(4) 




6.0 


1 


0.1652(5) 


0.14617(35) 


0.1243(7) 




5.0 


1 


0.1452(4) 


0.1283(5) 


0.1066(7) 


0.0955(5) 


4.5 


1 


0.1357(4) 


0.1187(4) 


0.1006(6) 


0.0895(10) 


4.0 


1 


0.1266(4) 


0.1095(5) 


0.0921(7) 


0.0813(4) 


3.5 


1 


0.1150(4) 


0.1019(4) 


0.0835(6) 





final uncertainty, we arrive at 

6.5 < gl < 12 (11) 
for the location of the supposed IRFP. 
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IV. MASS ANOMALOUS DIMENSION 

We derive the mass anomalous dimension from the 
scaling with L of the pseudoscalar renormalization factor 
Zp [aliUlllli- The latter is calculated by taking the 
ratio 



fp{L/2)- 



(12) 



fp is the propagator from the a wall source at the t = 
boundary to a point pseudoscalar operator at time L/2. 
The normalization of the wall source is removed by the 
/i factor, which is a boundary-to-boundary correlator. 
The constant c, which is an arbitrary normalization, is 
I/V2 in our convention. 

We present the calculated values Zp in Table III] and 
plot them in Fig. [51 Again , the slow running suggests 
that we may attempt to extract 7^ from the approximate 
scaling formula 



Zp{L)^Zp{Lo) ( 



(13) 



that is, from the slopes of the lines drawn in Fig.[6l These 
linear fits are analogous to Eq. (|5]): 



logZpiL/a) = Co + ci log 



L 



(14) 



The results are shown in Fig. [71 Similarly, we have also 
applied the linear fit with the smallest volume L = 6a 
removed, and we considered a fit function analogous to 

Eq. m, 



logZp{L/a) ^co + cilogL/8a + C2{\ogL/8af. (15) 



FIG. 6: The pseudoscalar renormalization constant Zp 
vs. L/a (Table |ll|). The crosses are from simulations with 
Pio — 0, P = 10.0 to 8.5. The circles are from simulations 
with /3io = 1: top to bottom, /? = 7.0 to 3.5. The straight 
lines are fits to each set of points at given (/3, /3io); the slope 
gives jm- 
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FIG. 7: Mass anomalous dimension 7((;^) from the linear fits 
shown in Fig. [S] plotted against g^{L = 8a). The squares are 
from the Pio = data while the circles are from /3io = 1. The 
line is the one-loop result. 
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FIG. 8: Comparison of different fit types. Empty symbols as in Fig. [7] Left panel: full symbols derive from linear fits in which 
the L = Qa points are omitted. Right panel: full symbols derive from fits to Eq. (|I5|1 . in which a log^ term has been added. 
Results are plotted against g^[L = 8a). The filled symbols have been slightly displaced to the right. 



In all cases the mass anomalous dimension at {L = 8a) 
is given by — ci. We show a comparison of the different 
fit types in Fig. [8l plotted against the running coupling 
g^{L = 8a). It is apparent that the result for ^mig^) 
is quite robust under variations in the fitting procedure. 
While some of the linear fits to all volumes give high x^, 
dropping the L — Qa points brings under control.^ 

A comparison of /3io = to /3io = 1 shows that there 
is some disagreement. Of the points obtained with /3io = 
1, only the weakest-coupling point is in agreement with 
the results of /3io = simulations. The two strongest- 
coupling points obtained with /?io ~ lie far above the 
line connecting the /3io — 1 points. The former originate 
from simulations near the strong-coupling transition of 
the /3io = theory. This is a lattice artifact, pushed 
off to stronger couplings by the introduction of /3io > 0. 
Thus the disagreement between the two sets of points 
should be settled in favor of the /3io = 1 points. 

All our fits show that 7^ departs from the one-loop 
line and levels off at 7„i ~ 0.45 at strong coupling. The 
highest point in Fig. |S]gives us the bound 7^ < 0.51. 



^ There are two /3io = 1 points with high where dropping L = 
6a leaves no degrees of freedom to the fit. One is the weakest- 
coupling point; it sits near the one-loop curve and agrees with 
the corresponding /3io = point, so we do not concern ourselves 
with it further. The other is the point at the strongest coupling, 
where x'^/dof =5/1 before dropping L = 6a. 



COMPARISON OF SU(2), SU(3), AND SU(4) 
THEORIES 



The present paper describes one of a set of three related 
theories, which differ only in their color content. The 
obvious way to compare these systems uses the language 
of large N . The 't Hooft coupling is A = g^N and in 
large N we expect to see collapse of data for different 
values of N onto a common function of A, up to 0{1/N) 
corrections. 

For theories with Nf — 2 fermions in the two-index 
symmetric representation, the renormalization group 
equation takes the form 



(16) 



d log ^ N 

where each of the terms on the right hand side is 0{N'~'), 



h 

N 



1 



167r2 
1 



167r2 



37V 



(17) 
(18) 



The large- limit of these systems has no IRFP in two- 
loop order, since 62 > 0. The one-loop mass anomalous 
dimension at A?^ = cxd is 



In 



6 



167r2 



A. 



(19) 



We collect our data from A^ = 2, 3, and 4 and present it 
in terms of A, taking g^ to be the Schrodinger functional 
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O A'=2 
O A'=3 
□ ^=4 

N=oo 




0.2 



FIG. 9: for the SU(2), SU(3), and SU(4) theories dis- 

played as a function of A = g'^N where is the Schrodinger 
functional coupling. The line is the lowest order large- A'' pre- 
diction, Eq. (I19|) . For SU(3), we have dropped two strong- 
coupling points calculated with = 0, since they are super- 
seded by Pa = 0.5 data; likewise for SU(4), where two /3io = 
points have been dropped (cf. Fig. [7]). 



coupling. For clarity, we use only the results of the linear 
fits, Eqs. dHl) and p^ . that do not drop L = 6a. 

We begin this time with 7™, plotted in Fig. [9l We 
have already noted that all three theories give results for 
7m that follow the one-loop line in weak coupling. Since 
the rescaled Casimir operators C2{R)/N for the three 
theories are similar (1, 10/9, 9/8) and close to the large- 
N value [1 -I- 0{1/N)]^ the ascending parts of the curves 
indeed almost coincide. What is remarkable is that all 
three theories fall off the one-loop line and level off in 
the neighborhood of 7™ = 0.4. This is a new regime of 
large- iV scaling behavior. 

The beta function, displayed in Fig. [TOl shows similar 
trends. In analogy with Eq. (O, we define 

and we plot it against u — 1/A. For our small values of 
iV, the leading correction to 61 is large, and so the weak 
coupling limits of the beta functions for different values 
of N do not coincide. We have shown the limiting behav- 
ior for each N on the right edge of Fig. [TUl along with the 
limiting value — 7/(247r^) [Eq. ([IT])]. Our results, along 
with their (possible) fixed points, march leftwards as N 



FIG. 10: Beta functions for the SU(2), SU(3), and SU(4) the- 
ories displayed a la large A'^ [Eq. (|20|) ]. The dashed lines at the 
right edge of the picture mark the weak-coupling limit as pre- 
dicted by the one-loop coefficient, Eq. p7p . The points with 
horizontal error bars mark each beta function's estimated zero 
(if it exists). 

increases. The iV = 2 IRFP lies at 1/A = 0.100(35). 
The existence of transitions for = 3 and 4 is more 
uncertain, but our fits put them at 1/A in the range 
0.044-0.067 and 0.021-0.038, respectively. It is certainly 
plausible to imagine that all three theories have IRFP's, 
and their location moves to ever stronger coupling as 
N increases. The two-loop beta function predicts that 
the IRFP has to disappear for N sufficiently large (at 
N ~ 37). We cannot rule out that by = 3 this has 
already happened. 

VI. CONCLUSIONS 

We began our study of these related models hoping to 
answer two questions: 

1. Do the systems exhibit walking, or do they possess 
an IRFP? 

2. How large does the mass anomalous dimension get? 

We were able to answer the first question cleanly for N = 
2: there is an IRFP. For larger N ^ the situation remains 
less clear: if there is an IRFP, it moves to ever stronger 
coupling as N increases, where simulations become ever 
more difficult. 

Qualitatively, the dynamics of all three models are 
dominated by their slow running. (This could have been 
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anticipated simply by examining the two-loop beta func- 
tion.) This makes the models quite different from con- 
ventional QCD with small Nf. Slow running was the key 
that allowed us to measure the mass anomalous dimen- 
sion. In all cases, 7™ remained less than about 0.45 over 
the observed range. This observation alone, independent 
of the existence of an IRFP, renders the theories unsuit- 
able as candidates for phenomenologically viable walking 
technicolor models. 

We are struck by the fact that, for all three theories, 
the two-loop beta function, expressed in terms of the 
Schrodinger functional coupling, does an excellent job of 
reproducing our observed beta functions. At the same 
time, the behavior of the mass anomalous dimension is 
quite different from the perturbative prediction, and is 
universal — for all three theories 7™ falls off the one loop 
formula and becomes (nearly) constant, independent of 
the value of the renormalized gauge coupling. 
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Appendix A: nHYP smearing for SU(4) 

1. Smearing and normalization 

The formulas of nHYP smearing, introduced in 
Ref. [13], are given in detail in Ref. [s^ for the SU(2) 
and SU(3) gauge groups. The smeared links Vnfj. are 
constructed from the bare links t/„^ in three consecutive 
smearing steps via intermediate fields V and V according 
to 



(Ala) 



(Alb) 



(Ale) 

I 

matrix V is given by 

V = n{n''n)-^/^. (A2) 

This requires calculation of the inverse square root oiQ = 
il^il, which is a positive Hcrmitian matrix. Presuming Q 
to be non-singular, the Cayley-Hamilton theorem allows 
us to write Q'^^"^ as a polynomial in Q, 

Q'^'^ = fo + fiQ + hQ^ + hQ^- (A3) 

We use a Jacobi algorithm to solve for the eigenvalues gi 
of Q. (We will not need derivatives of gi so they need 
not be found analytically.) By writing Eq. (|A3[) in the 



±p^!^,p 

(l-a3)C/nM + Y Yl U^^vUn+f,,pUl^^^^ . 
±ri^p,iy,p 
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eigenbasis of Q we obtain a linear system for fi, 

-1/2 \ 



1 90 90 9o 
1 91 9l 9i 
1 92 gl 92 







] 













9q 

-1/2 

9i 

-1/2 
92 



(A4) 



This system can also be solved numerically; since, how- 
ever, we will need to differentiate the result, we solve it 
analytically by inverting the Vandermonde matrix. The 
solutions fi are rational functions of = or, more 
conveniently, of the symmetric polynomials 



u = ro + ri + r2 + 
V = rori + rQr2 + Tors + rir2 + rirs 
+r2r3 

w = rorir2 + rorir^ + ror2r3 + rir2r3 
X = rorir2r3. 

Denoting a common denominator by 



we have fi 
No 

^2 ^ 

N3 : 



A = a; (^u^x — uvw + vo^ 
Ni/A, where 
v?wx — u (— u^a; + vuP' + 

u^x - 2u'^vw + u{v^ + 2w'^) 
—w (w^ + x) 

—u^v + v^w — u i^x — 2u^) — 



(A5a) 

(A5b) 
(A5c) 
(A5d) 

(A6) 



X — vwx 



2vw 



(A7a) 

(A7b) 
(A7c) 
(A7d) 



This completes the calculation of the quantities needed 
to normalize O. 



2. Force in molecular dynamics 



We follow still the derivation in Sec. 3 of Ref. 134 



which is based on Morningstar and Peardon [42 1. 

The force is the derivative of the effective action with 
respect to simulation time r. The first step is to note 
that the fermionic part of the action includes only the 
fat links Ki/x, so 

A5^g^Retr^^^Retr(E„^K^). (A8) 

One proceeds to apply the chain rule repeatedly to 
via Eqs. (|Aip until one arrives at derivatives t/„^ of 
the thin links.'' The only factor in the chain rule that 



In fact the chain rule is first applied to the change of representa- 
tion from decuplet to fundamental; this is followed by the fat-link 
chain rule. 



depends on the group comes from the U(7V) projec- 
tion (jA2p . which appears at every level of smearing in 
Eqs. (|Al|) . In order to write a derivative V in terms of (l, 
we use the Cayley-Hamilton formula (jA3l) [cf. Eq. (3.10) 
of Ref. [Hi, 



RetrSF = Rctr 



Rctr (g^^/^sii) + tr (sr2)/o + tr (Qsrj)/i 

+tr (Q2sri)/2 + tr {Q^J:n)h + fi tr (S^iQ) 

+f2 tr [{^nq + QEn)Q] + h tr [(srjg^ 



Upon defining the traces 



n + l 



n+l 



one can write 



71=0 

where bin = dfi/dcn- Then 

RetrSy = RetT {Q-^^^J:n) + RetT AQ, 
where 



(A9) 



(AlO) 



(All) 



(A12) 



A = ^tr(B„Er2)Q" + /ii;r2 + /2(i]f70 + oi]f7) 



(A13) 



with Bn = bon + binQ + b2nQ^ + hnQ^ ■ Now wc differ- 
entiate g = to obtain finally 



RetrSF = Rctr 



(Q-^^^T, + An'' + A^n'')n] . (ai4) 



The derivatives bij are calculated via the eigenvalues 
gk through the chain rule. 



l^.. ^ ^ ^fi ^9k 

9cj ^ dgk dcj ' 



(A15) 



The matrix dgk/dcj is the inverse of the Vandermonde 
matrix dck/dgj = (gj)'^, so that we still don't need an 
analytical expression for the eigenvalues. The derivatives 
dfi/dgk can be calculated directly from the above ex- 
pressions for fi. The final result for bij can, like the fi, 
be written as rational functions of the symmetric poly- 
nomials u,v,w,x. The expressions are lengthy and not 
particularly illuminating '43] . 
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Appendix B: Tests of soft gauge actions 

In lowest order, the action of Eq. ([3]) is just a quadratic 
form in the vector potentials of the thin and fat links. 
In momentum space the fat link's gauge field B^{q) is 
related to the thin link's gauge field A^{q) through a 
form factor h^i,{q) whose specific form depends on the 
particular definition of the fat link, 

B^{q)^Y.~^^M)AM- (Bl) 

In general, - 1 + 0{a'^q^) + 0{a^q^) + ■■■. 

Let us generalize to N colors, with an ordinary pla- 
quette term made of thin link plus a closed loop made of 
fat links in representation R of the gauge group. We can 
write the quadratic gauge action as 

^ffo Jpp' 

X [{l-s)A,{p')D'^^p)AM 
+sB^{p')D';^{p)B,ip)] , (B2) 

where is the thin link kernel and D^'^^ is the fat link 

kernel. (We have suppressed color indices.) In practice, 
we will use a plaquette fat link term, so D^^^ — _D"^™. 
The bare coupling is 



\ = ^T{f) + ^T{R), 
9o N dn 



(B3) 



where T(/) = 1/2 is the trace normalization of the fun- 
damental representation while T{R) is that of the repre- 
sentation R; d{R) is the dimension of the latter. We also 
define s = x/{l -\- x) where 



(3RT{R)/di, 

mf)/N 



(B4) 



In SU(4) with decuplet fat links, x = ^i3r/(3. 

Upon using Eq. (|Bip . wc find the free gauge boson 
action to be 

1 



So=-7r2 I {27rr5\p + p')[A^{p')D,,{p)A^{p)], 

(B5) 

where 



D,, = (1 - s)Dl';r 

After adding a gauge fixing term 
1 ^ f 1,- 



Sg[ — 7^~2 / ^ / Tr— fc^/Cj/A^( fc)yli/(fc), 

^■9o „ 



(B6) 



(B7) 



where A:^ — (2/a) sin(aA;^/2), the gauge boson propaga- 
tor is found by inverting the field equation. 



E 



1: 



(B8) 



The case of a Wilson action and APE smearing can 
be treated analytically. In that case both the action and 
smearing are built out of projectors. 



= S 



fc2 ' 



/If 



fc2 



The Wilson gauge action is 

Dfivik) — k P^jy, 
and the APE smearing term is [i^l 



(B9) 



(BIO) 



(Bll) 



with f{q) = 1 — {a/6)q'^; a is the smearing parameter. 
This means that the soft APE-smeared Wilson action is 



D^,{k) = (i-s + sf)k'p;^, 

and thus the propagator is 



Gf_iu{k) — j- 



1 



1-S + sf 



-p' 

2 fJ-f 



(B12) 



(B13) 



More complicated actions and smearing do not have this 
projector form, and the final result for G^u{k) is not il- 
luminating. 

So much for the gluon propagator. The second ingredi- 
ent needed for perturbative calculations is the fcrmion- 
gauge boson vertex. In a fat action with unitary links 
it is simply the smearing form factor multiplied by the 
unsmeared vertex F^, 



(B14) 



All Feynman graphs with an internal gluon line are built 
of the combination of terms 

Gpu{k) — Y ^iG p^{k)T I, — Fy [K^G^pKp\T^ (B15) 

and so for all practical purposes, the gluon propagator 
Gfj,,y{k) is just replaced by the combination h^^G^ph^p. 
For APE smearing on top of the Wilson gauge action, 
this replacement can be performed analytically, and the 
transverse part of the propagator times form factors is 
just the usual Wilson propagator, rescaled by the factor 



l-s + sp{qy 



(B16) 



Let us recall some features of perturbation theory with 
fat links. Consider first the usual smeared fermion action 
with unsmeared gauge action, s = 0. The good features 
of this action arise because the form factor in Eq. (|B16P 
suppresses the large-g part of loop integrals. This shrinks 
the size of tadpole graphs and lowers the size of one loop 
matching factors. A positive weight s increases the form 
factor at large g, and a negative s effectively amounts to 
increased smearing. 
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FIG. 11: Additive mass renormalization [see Eq. HB17|l ] for 
nHYP clover fermions with csw = 1 as a function of the 
relative weight s of thin- and fat- link terms. 



The case s = 1 is peculiar in that the effects of smear- 
ing are completely undone. What is happening in this 
case is that the transformation A B is just a field re- 
definition: we are doing Monte Carlo with the fat link 
everywhere as the fundamental gauge field. Only if the 
gauge action is smeared differently from the vertex can 
smearing achieve its goal. 

One way to reduce the value of the effective gauge bo- 
son propagator at large q is to make s negative. Then 
the fat term amounts to a further smearing of the effec- 
tive quark-gluon vertex. This might be a feature worth 
exploring in low-TVj situations. 

As described in the main body of the paper, we have 
chosen s > for a nonperturbative reason, and so we 
should check if this choice is benign from the point of 
view of perturbation theory. We show the results for two 
observables to check this point. 

Fig. [TT] shows the additive mass renormalization for 
nHYP clover fermions with csw = Ij parametrized as 



5m = ^2T{R)So. 

477 



(B17) 



(The choice of factors is taken to make the case of funda- 
mental representation fermions most transparent.) The 
s = 1 case reduces to the thin link value. The undoing 
of smearing is quite abrupt, and nearly any value away 
from s = 1 produces a large suppression of Sm. 

In perturbation theory, the tree level static potential 
between point sources is Coulombic, up to lattice arti- 
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FIG. 12: 4:TTrV{r), the tree-level static potential due to a soft 
Wilson action and seen by nHYP-smeared sources, at four s 
values: (a) s = 0.7, (b) s = 0.2, (c) s = 0, (d) s = -0.5. 



facts: 



(B18) 



When the sources are themselves nHYP-smeared, the 
short distance part of the potential is softened. This 
is a usual consequence of smearing . Fig. [TH shows 
A7TrV{r) for soft actions with several choices of s. It 
appears that from the point of view of perturbation the- 
ory, the soft gauge action does not obviously harm useful 
properties of the lattice action. 

Finally, the simulations we are performing with the 
soft action are done at rather strong coupling. In this 
regime, perturbation theory does not reliably predict the 
observed additive mass renormalization. This is shown 
indirectly in Fig. [131 a plot of Kc against the bare cou- 
pling 8/(7q = /3 + ^Pio- Different plotting symbols show 
different values of fiiQ. In weak coupling the observed 
additive mass renormalization qualitatively agrees with 
what is shown in Fig. 1111 larger /3iq / (3 produces larger 
additive mass renormalization. However, in strong cou- 
pling the effect reverses. One cannot help speculating 
that there is a connection between a large value of Kc 
and a nearby first order transition, and that when Kc 
falls, the transition has moved away. 



Appendix C: Ensembles 

Our algorithm is the hybrid Monte Carlo (HMC) al- 
gorithm of Duane and Kogut [46l - l48j . We accelerate the 
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FIG. 13: Dependence of Kc on bare coupling S/go = /3 + ^/3io, 
for various values of /3io. 



molecular dynamics integration with an additional heavy 



pseudo-fermion field as suggested by Hasenbusch 
with multiple time scales [5Q] . and with a second-order 
Omelyan integrator [5l[. Lattice sizes range from 6^ to 
16^ sites. 

We list in Table Hill the values of (/3, /3io, 'tc) and the 
number of trajectories run at each volume, along with 
the length of the trajectories and the acceptance. Poor 
acceptance forced us to shorten the trajectory length in 
many cases from the usual value of 1. 

The observables we measure are the (inverse) 
Schrodinger- functional running coupling, l/g^, and the 
pseudoscalar renormalization factor, Zp. (We measure 
Zp on the same configurations used to determine l/ff^-) 
Both of them turn out to have long autocorrelations. We 
monitored and controlled this problem by running 4 or 
8 streams in parallel at each (/3,/3io) and L. After an- 
alyzing each stream separately, we fit the results of the 
streams together to a constant. We demanded that the 
X^/dof of the constant fit not exceed 6/3 for 4 streams, 
or 10/7 for 8 streams. For the largest volume L = 16a 
at (/3,/3io) — (5,1), we were not able to overcome the 
autocorrelations in 1/g'^. Wc therefore omit this point 
from the analysis of the running coupling. The autocor- 
relations in Zp, on the other hand, did allow a consistent 
determination, and thus we keep this point in the analysis 
of the mass anomalous dimension. 
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TABLE III: Ensembles generated at the bare couplings 
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